NPS ARCHIVE 
1963 

JENISTA, 3 . 



APPLICATION OF THEORETICAL DESIGN METHODS 
TO THE STUDY OF PERFORMANCE LIMITS ' 

OF AIRFOILS IN CASCADE 

JOHN E. JENISTA 




library 

NAVAL PCSTCRADUATE SCHOOL 
MONTEREY. CALIFORNIA 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



APPLICATION OF THEORETICAL DESIGN METHODS 
TO THE STUDY OF PERFORMANCE LIMITS 
OF AIRFOILS IN CASCADE 



# 3$C # # # 



John E. Jenista 



APPLICATION OF THEORETICAL DESIGN METHODS 
TO THE STUDY OF PERFORMANCE LIMITS 
OF AIRFOILS IN CASCADE 



by 



John E. Jenista 

// 

Lieutenant Commander, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

AERONAUTICAL ENGINEERING 

United States Naval Postgraduate School 
Monterey, California 



1963 



'A"\ 








































LIBRARY 

U S. NAVAL POSTGRADUATE SCHOOL 
r MONTEREY. CALIFORNIA 



APPLICATION OF THEORETICAL DESIGN METHODS 
TO THE STUDY OF PERFORMANCE LIMITS 
OF AIRFOILS IN CASCADE 
by 

John E. Jenista 

This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

AERONAUTICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



A method of designing airfoils in cascade by means of 
conformal transformations is discussed. Design of these air- 
foils is regulated by five independent input parameters, with 
solutions obtained by digital computer. A large number of 
cascades generated by this method are compared. To evaluate 
the limits of performance, a parameter to indicate the ten- 
dency toward flow separation is introduced, with a limiting 
value established and verified. Proper solidity is shown to 
be of great importance in achieving low values of this separ- 
ation parameter. The value of proper solidity for a given 
blade thickness is shown to be relatively independent of 
fairing shape. To increase performance, reducing blade thick- 
ness with a corresponding increase in solidity is shown to 
yield much greater improvement than changes in fairing shape. 
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1.0 Introduction 

This study is an extension of a line of investigation 
originally conducted by Professor Theodore H. Gawain of the 
U. S. Naval Postgraduate School. The original investigation, 
sponsored by the Aerojet General Corporation, was specifically 
aimed at increasing the pressure rise per stage of axial flow 
pumps. The intended application was to reduce the size of 
propellant pumps for large liquid fuel rockets. The require- 
ments were for maximum pressure rise per stage, regardless 
of efficiency. This, in turn, demanded positive and accurate 
control of the pressure distribution around each blade, with 
particular emphasis on the adverse pressure gradient experi- 
enced by the boundary layer. 

Flow of a real fluid through an axial flow compressor or 
pump involves some formidable problems. Besides the effects 
of viscosity and turbulence; tip losses, centrifugal force 
fields, and other effects complicate the problem considerably. 
As is usually the case in engineering however, a simplifica- 
tion of the problem may yield useful results. If the stream 
surfaces of the flow are assumed to be concentric cylinders, 
then a cylindrical cross section of the flow may be repre- 
sented by an infinite linear cascade with two-dimensional flow. 
If all velocities are considered relative to the blade row, 
then the solution for the rotor and the stator have identical 
mathematical form. Following the usual aerodynamic practice, 
viscosity may be neglected outside the boundary layer. This 
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leaves only turbulence as a factor preventing the use of poten- 
tial flow considerations. In order to use the powerful analy- 
tic methods of potential flow, turbulence has been neglected. 
This does not mean however, that the data presented is invalid. 
Actually, the potential flow assumption is used to obtain the 
pressure distribution only. The relation between pressure 
distribution and separation was obtained from experimental 
data from actual turbomachines, and must therefore include the 
effects of turbulence. 

From the foregoing considerations, it may be seen that 
the study of cascade design has application to the whole field 
of axial flow turbomachinery. The problem of cascade design 
is much more difficult than the problem of isolated airfoil 
design. In the early days of airfoil design, much experimen- 
tal data was obtained on a bewildering variety of shapes. 
However, there was no reliable method of predicting the per- 
formance of an airfoil before it was tested. It was only 
after the NACA (now NASA) discovered how airfoil performance 
depended on shape, specifically on thickness, camber and fair- 
ing distribution, that any real progress was made in the field 
of airfoil design. 

The present state of the field of cascade design is that 
there is a large volume of experimental data available on a 
bewildering variety of cascades. Many efforts have been made 
to extrapolate the theory of isolated airfoils so that it will 
apply to cascades. However, at present there is no reliable 
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theory which covers the field of cascade performance in the 
way that NASA's theory covers the field of airfoil performance. 
The parameters of thickness, camber, and fairing distribution 
are still involved, plus a few more. The relations and in- 
teractions between these parameters are so complex, that it 
is very difficult to consider any aspect of cascade design 
separately. 

When considering cascades as complete entities, the inter 
action effects are all present, and need not be corrected for. 
However, experimental efforts in this direction have been nec- 
essarily limited in scope, due to the effort and expense in- 
volved in constructing and testing large numbers of cascades. 
Theoretical methods, such as the one outlined in this investi- 
gation, provide the means of considering complete cascades 
over a wide range of design parameters. This study was under- 
taken from the standpoint of discovering areas of possible 
improvement in compressor performance. 

Most of the efforts of the NASA in the field of compre- 
sor design have been for maximum efficiency with moderate 
output. In some cases, it may be required to attain maximum 
output, with perhaps moderate efficiency. The gains to be 
expected from these considerations are not of large order how- 
ever, as a cascade develops its best efficiency not far from 
the point of stalling the blades. 

The basic problem, then, is composed of two parts. First 
it is necessary to devise a method for generating blade shapes 
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with the corresponding pressure distributions. For this pur- 
pose, the powerful methods of complex variables and conformal 
transformations are utilized. The method developed by Pro- 
fessor Gawain generates infinite cascades of airfoil shapes, 
based on five arbitrary parameters. Generally, the transfor- 
mations are similar to the Joukowski transformations with 
additional requirements. 

The second part of the problem is to establish a correla- 
tion between the pressure distribution around each blade and 
the occurence of separation. This problem must, of course, 
depend heavily on experimental data. Fortunately, work has 
been done in this area and some practical working limits have 
been established. See Reference 1. Based on this data, the 
simple parameter ( V max /V 2 )^ has been chosen as an indication 
of the tendency toward separation. 

The author wishes to express his sincere appreciation to 
Professor Theodore H. Gawain of the U. S. Naval Postgraduate 
School for his assistance, guidance, and encouragement during 
the course of this investigation. 
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2.0 The Method of Cascade Design 



2,1 introduction 

Methods of designing airfoils by purely mathematical means 
have been known since the early days of aviation. Joukowski 
was the most notable pioneer in this field. He, and others, 
noted the possibility of designing cascades of airfoils by 
similar methods. Since the early applications of cascade de- 
sign were to multi-winged airplanes, not much work was done in 
this field. As pointed out earlier, the application to turbo- 
machinery makes this field of considerable interest today. 

In considering the flow through cascades, the shape and 
position of any individual blade will affect the flow around 
that blade. In addition, the shape and position of all the 
other blades will also affect the flow around that same blade 
at the same time. Therefore, more parameters must be used in 
cascade design than are used in isolated airfoil design. In 
addition to blade shape, the relation of the blades to each 
other, and their relation to the flow must all be described. 

The method developed by Professor Gawain uses five independent 
parameters to describe an individual cascade. 

The solution for blade shape and velocity distribution 
corresponding to chosen values of the five parameters is found 
by using relations of potential flow and conformal mapping in 
three separate planes. These have been termed the near circle 
plane, the circle plane, and the cascade plane respectively. 
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In reading through the following sections, which describe 
the general method of generating cascades, reference may be 
made to Figure 5> which illustrates a typical contour in each 
of the planes referred to above. A summary of the mathematical 
development of the method is given in Appendix A. A more de- 
tailed analysis may be found in Reference 2. 

The solution of these equations, many of which are trans- 
cendental, is an impossible task by hand calculation. The 
entire problem has been arranged in FORTRAN language for solu- 
tion by a digital computer. The CDC 1604 computer at the U. S. 
Naval Postgraduate School was used for this investigation. 

The basic computer program, along with later modifications is 
included in Appendix B. 

The following sections describe the general method of 
generating cascades from five input parameters. Each parameter 
is described in detail as it enters into the calculations. 

The parameters are: 

ft = Mean flow direction 

B 0 - Turning parameter 

O'' - Equivalent flat plate solidity 

S' 85 Thickness parameter 

/< = Shape or symmetry parameter 

2.2 The Flow Parameters 

Flow through a cascade is usually described in terms of 
two velocities. The first velocity is that far enough up- 
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stream of the cascade so that local effects due to the cascade 
are not evident. This is called the inlet velocity, and the 
angle it makes with the reference axis is called Similarly, 

the velocity far enough downstream so that local effects have 
damped out is called the outlet velocity, and its angle with 
the reference axis is called The magnitudes of these two 

velocities are not independent, since the laws of continuity 
require that the component perpendicular to the cascade axis 
be equal for the two cases. This concept allows velocities to 
be conveniently nondimensionalized, so that only angles are 
important . 

Figure 1 is an illustration of two typical cascades. An 
appropriate velocity vector diagram is shown for each. Note 
that only two parameters are necessary to describe the complete 
diagram. Instead of the angles and /S 3 it was more conven- 
ient to use two different parameters. Consider the following 
expressions: 

Tan (3 * £(Tan /3^ + Tan (3 2 ) 

B c = (Tan fa -Tan /S 2 ) 

The first expression defines the angle /?, or mean flow direc- 
tion, which is the first input parameter. The second expres- 
sion defines the turning parameter B c , which is the second 
input parameter. The relationship between all these quantities 
is shown in the diagram on Figure 1. 

The starting point of the calculations is the near circle 
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plane. In this plane, the complex potential for a series of 
doublets spaced at intervals of 2 Tt along the y axis is com- 
bined with the potential for uniform flow at the angle ft with 
the x axis. The velocity of this flow is set at unity, which 
then becomes the reference for all the nondimensional veloci- 
ties derived later. If only one doublet were used, the resul- 
ting potential would describe the flow without circulation 
around a single circular cylinder. In this case, with an in- 
finite series of doublets along the y axis, the resulting 
potential describes the flow without circulation around an 
infinite series of nearly circular bodies spaced along the y 
axis. These bodies are not actually circles because the pre- 
sence of the infinite stack of doublets distorts the flow 
around any individual doublet, such that the flow contour is 
a slightly flattened circle. The amount of distortion is re- 
lated to the strength of the doublets and the spacing between 
them. In this case, the spacing between doublets is fixed at 
2 Tf , so doublet strength is what determines the relation 
between spacing and size of the bodies in the near circle plane. 

2.3 The Solidity Parameter 

The strength of the doublets is fixed by the input para- 
meter CT , or the equivalent flat plate solidity. This para- 
meter is closely related to, but not identical with, the ordinary 
solidity c/S, or the chord-pitch ratio of the blade row. The 
distinction can best be illustrated by an example. Generally, 



any row of arbitrary blades of solidity c/S may be transformed 
conformally into a related row of flat plate vanes. The chord- 
pitch ratio of the latter is then denoted by O' . It follows 
that if two arbitrary cascades have the same value of (J* , 
then either may be conformally transformed into the other. But 
if the value of CT is different for each, then no such trans- 
formation is possible. (T" then, may be regarded as a kind 
of generalized solidity, of perhaps deeper aerodynamic signi- 
ficance than simple geometric solidity. 

The relation between values of 0"~ and geometric solidity 
is directly influenced by blade thickness. Very thin blades 
have geometric solidities almost exactly equal to CT . With 
thicker blades, the geometric solidity is always less. The 
reason for this may be illustrated by examples of two extremes. 
If a contour in the near circle plane is transformed confor- 
mally into a flat plate, the chord length of the flat plate 
will be nearly twice the diameter of the original contour. 
However, if no transformation is made at all, an "airfoil" of 
100$ thickness results, whose chord length is equal to the 
diameter. Both cascades would have the same value of (T* , but 
the geometric solidity of the latter case would be half that 
of the former. Between these two extremes, finite thickness 
airfoils in cascade have solidities less than . 

2.4 The Turning Parameter 

The turning parameter, B c , is a measure of the change in 
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direction of the flow passing through the blade row. It is 
actually a nondimensional form of the downwash velocity induced 
by the blades. There is a very simple relation between B 0 and 
the circulation around each blade. In fact, for a constant 
mean flow angle, this relation is a simple direct proportion. 
This convenient state of affairs is the reason why the two 
flow parameters were expressed in their present form. 

In isolated airfoil theory, the method of introducing lift 
or circulation on an airfoil, is to add to the complex poten- 
tial, a term representing a vortex of suitable strength. How- 
ever, mere addition of a vortex will produce a distortion in 
the contours in the near circle plane, since these contours 
are not exactly circles. In addition, an infinite stack of 
vortices is necessary to produce circulation around each of the 
blades in the cascade. To avoid these resulting undesirable 
perturbations in what will ultimately become the airfoil shape, 
and indirect method of adding circulation is used. 

In this method, a transformation is made from the near 
circle plane to the circle plane. An equation is developed 
which transforms the flow in the near circle plane into the 
flow around a perfect circle, centered on the origin. The cir- 
culation term may then be introduced without fear of distor- 
tion. The transformation is then made back to the near circle 
plane. Since the same equation is used to transform both ways, 
the original flow contour remains unchanged; however, the 
velocities at each point are now different, due to the circu- 
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2.5 The Blade Shape Parameters 

The final transformation is from the near circle plane 
to the cascade plane. This is done by the Joukowski method 
of shifting the axes slightly in such a way that the flow con- 
tour transforms into an airfoil shape. Two arbitrary para- 
meters were chosen to control the blade shape. One of these 
is called the thickness parameter, and is given the symbol <S . 
This quantity is generally related to the thickness, and also 
to nose curvature, but not in any fixed relationship. Large 
values of £, produce thick airfoils, and small values produce 
thin airfoils. Figure 2 is an illustration of several airfoils, 
showing the effect of < £> . The limiting value of 0 for £ does 
not however produce an airfoil of zero thickness. Due to the 
slight distortion of the contours in the near circle plane, 
6=0 may produce the mathematically significant, but physi- 
cally useless case where the "lower" surface of the airfoil 
crosses through the "upper” surface. 6 then, can be small but 
not zero, and the airfoils produced can be thin but not of zero 
thickness. 

The other parameter is the shape or symmetry parameter K. 
The effect of variations in K may be seen in Figure 3. Gener- 
ally, K“0 produces a blunt leading edge, and cusped trailing 
edge similar to a Joukowski airfoil. K=1 will produce a sym- 
metrical shape which is generally a cambered ellipse. As K 
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increases from zero to one, the point of maximum thickness 
moves from the quarter chord to the half chord point, and the 
trailing edge rounds off* Values of K larger than one will 
produce airfoils facing the wrong way (Trailing edge into the 
flow) and hence must not be use do 

2 .6 Stagnation Points 

During the transformations made in this analysis, the 
stagnation points serve as reference points to locate the re- 
lation of the contours in the various planes* For this reason 
there is always a stagnation point at the leading and the 
trailing edge* Having a stagnation point at the leading edge 
is sometimes called the Theodorsen condition. This means that 
in all the cascades developed by this method, the airfoils are 
operating at the ideal angle of attack* Variations in lift 
then, are achieved by variations in the camber of the blades* 
Figure 1 graphically illustrates this fact* Two typical 
cascades are drawn, with the principal variation being in turn- 
ing parameter, B c * The values of £, were adjusted to give 
the same blade thickness* The difference in camber may be 
clearly seen* 

The condition of a stagnation point at the trailing edge 
is called the Kutta condition* This fact has particular sig- 
nificance when considering the blade shapes where K=l* As 
pointed out previously, these shapes are generally cambered 
ellipses* It is a well known fact that flow around an ellipse 
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of a given thickness ratio produces a lower velocity peak than 
the flow around any other shape of the same thickness ratio. 
This class of blades then should exhibit the lowest loading 
attainable for any possible fairing shape. Of course, any 
attempt to duplicate this fairing shape in an actual cascade 
is doomed to failure because of the rounded trailing edge. 
However, this case still has significance in two respects. 
First, the performance of a cambered ellipse, with the Kutta 
condition mathematically imposed, represents a useful limiting 
ease of performance due to fairing shape. Secondly, a blade 
similar to this, with modifications to produce a sharp trail- 
ing edge, would have very nearly the same performance. The 
critical area of flow near the leading edge would not be af- 
fected much by modifications at the rear of the airfoil. 

Figure 4 is a plot of the pressure distribution around 
the airfoils of Figure 3» The above premise is verified, as 
the negative pressure peak can be seen to decrease as the 
value of K is increased. Note that the pressure distribution 
aft of the peak is essentially linear to the trailing edge. 
Reference 3, p* 354 states that this type of pressure distri- 
bution is typical of decelerating cascades. Note also that 
the slope of this line becomes less as K is increased, indi- 
cating a less adverse pressure gradient in the boundary layer. 
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3.0 Prediction of Separation 



3*1 Pressure Gradient 

The phenomenon which ultimately limits the attainable 
pressure rise of a cascade is separation. Flow separation 
occurs when flow in a boundary layer experiences an adverse 
(positive) pressure gradient. This adverse gradient tends to 
decrease the velocities in the boundary layer. Ultimately, 
the point is reached where the velocities close to the surface 
have been reduced to zero. At this point, flow separation 
occurs. Beyond this point, some velocities become negative, 
and a region of reverse flow exists. The resultant effects 
cause high drag and loss of circulation around each blade 
element. Many experiments have verified that any adverse 
pressure gradient, however small, will cause separation if it 
acts long enough. The concept of a gradient and the distance 
over which it acts, leads one naturally to suspect pressure 
difference as a significant parameter in predicting separation. 
The gradient, integrated over the distance is mathematically 
equal to the difference in pressure. It was previously shown, 
furthermore, that the pressure distribution over typical air- 
foils in this system is approximately linear aft of the nega- 
tive pressure peak. 

3.2 Velocity Considerations 

From another point of view, one might examine the velo- 
cities of flow through a cascade. From the inlet velocity 
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described previously, velocity is brought to zero at the 
front stagnation point, accelerated to a local maximum some- 
where on the surface of the blades, reduced toward zero at 
the trailing edge, and finally settles down to the outlet 
velocity far behind the cascade. Besides the inlet and outlet 
velocities, only one other velocity is really significant, and 
that is the local maximum, V max . From these three quantities, 
only two independent dimensionless parameters may be construc- 
ted. If the inlet velocity and V max were used, this parameter 
would be related to the flow around the forward part of the 
airfoil. Since this flow involves a negative pressure grad- 
ient, it is probably not very significant with respect to 
separation. Similarly, the quantities V max and the outlet 
velocity, called V 2 , are related to the flow around the rear 
part of the airfoil, which involves a positive pressure grad- 
ient. These two parameters then, combined in some dimension- 
less manner, should act as an indicator for the tendency 
toward separation. 

3 .3 Diffusion Factors 

From its many experimental studies of flow through cas- 
cades, the NASA has arrived at essentially the same conclusion. 
They have used several dimensionless combinations of V max and 
V 2 to define blade loading limits. One of these is the 
"Local Diffusion Factor” which is: 
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(V - V ) 

D-, = max 2 

V 

max 

This factor is described in Reference 1 and supporting data 
indicates that separation is likely to occur for values of 
greater than 0.5. 

Since the velocity distribution around the blades must 
be known to find values of , it is necessary to use experi- 
mental data, or an analytic al treatment such as the one des- 
cribed in this investigation. To get around this difficulty, 
the NASA also defines an approximation, D, which may be ob- 
tained from more easily determined quantities. D, is called 
the "Diffusion Factor" and is given by: 

D - 1- Cos + Cos (Tan ftl - Tan/92) 

cos p z 2 y s 

It is seen that this quantity can be determined from the inlet 
and outlet velocities, and the solidity of the blade row. The 
supporting data indicates that separation is likely to occur 
for values of D greater than 0.6. 

Since the methods used in this investigation furnish all 
of the flow quantities described above, as well as the com- 
plete velocity distribution around each blade, the opportunity 
to compare the two NASA parameters is presented. As will be 
described later, the correlation between these two parameters 
is not too good for the solidities considered here. 
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3«4 The Boundary Layer Loading Parameter 

From the above development? it may be seen that the para- 
meter to indicate separation should involve the pressure ratio 
across the rear part of the airfoil, and be some dimensionless 
combination of the quantities V 2 and V max <> The NASA has in- 
dicated that most compressors operate in the range where 
Reynolds number effects are insignificant <> Therefore, it is 
felt that no particular inaccuracies result from omitting 
consideration of Reynolds number 0 

In order to avoid using one parameter for pressure dis- 
tribution, and another for the tendency toward separation, 
the quantity S was selected where: 

s - ( W v 2 ) 2 

This quantity is referred to as the boundary layer loading 
parameter, and is quite versatile in application.. It is 
directly related to the pressure coefficient, and in fact only 
differs from that quantity by unity., S also has a constant 
functional relationship with the NASA parameter D^o The value 
S = 4.25 corresponds to the value of of 0,5 1 and has been 
selected as the limit loading in the boundary layer for blades 
in this investigation,, 
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4.0 Results and Discussion 

Theory and experiment both point toward increasing 
solidity as the method of obtaining increased pressure rise 
per stage. As solidity is increased, the lift coefficient on 
each blade is decreased, and hence the loading on the boundary 
layer is decreased also. However, as the solidity is increased, 
there is an adverse effect due to the blockage caused by the 
finite thickness of the blades plus the boundary layer. 

Finally, the point is reached where the beneficial effects of 
decreased blade loading are just counterbalanced by the ad- 
verse effects of blockage. The quantitive determination of 
this point was one of the major objectives of this investiga- 
tion. 



4.1 Organization and Scope 

Since actual cascades are designed on the basis of blade 
thickness and geometric solidity, it was felt that the results 
of this investigation ought to be organized on this basis, 
rather than on the basis of the mathematical parameters <£✓ and 
O'" . The value of blade thickness is not directly predictable 
from the input parameter £ . Therefore, the basic computer 
program was modified to use an iterative scheme to achieve the 
desired values of blade thickness. Details are given in 
Appendix B, which describes the basic computer program and the 
modifications. Seven blade thicknesses were selected to cover 
the field of possible blades. These ranged from the limiting 
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case of extremely thin blades to ^4% thickness blades in four 
percent increments., 

It was decided to conduct this investigation using only 
one value of average flow angle /3. This value was chosen to 
be 45°, which was a somewhat arbitrary decision, but not with- 
out some sound theoretical basis. It is stated in Reference 
4 that maximum efficiency may be obtained when the average 
flow angle is near 45°, and many compressors are built accord- 
ingly. 

In spite of reducing the parameters to four by the above 
decision, over 1470 separate cascades were generated and ana- 
lyzed for this investigation. Additional parameters could 
have been introduced to give more control over the blade shape, 
but they would have greatly magnified the problem. As the 
data will show, additional blade shape parameters would not 
have affected the significance of the results very much. 

4.2 Series 1 Plots 

The first plots of data from the computer run were made 
of boundary layer loading vs. solidity, for airfoils of con- 
stant blade thickness. Five values of the shape parameter K 
were used, and data was obtained at six values of turning 
parameter, B c . This meant that there were 30 graphs in series 
one. Figures 6 through 9 are representative graphs from this 
series. 

It was noted with satisfaction that the curves showed a 
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pronounced minimum loading as solidity was increased. The one 
exception was the case of the extremely thin airfoils, for 
which the curve appeared to be asymptotic to some low value of 
loading. This is not unreasonable, since the blockage effect 
would be very small for this case. The case of thin airfoils 
then, represent? a sort of limiting value for cascades. The 
point at which minimum loading occurred for a given thickness 
airfoil was called the optimum solidity. 

It may be noted that while the value for loading at an 
optimum point could be determined with accuracy, the corres- 
ponding solidity for that point was subject to some uncertainty. 
Therefore, another point was introduced where the solidity 
could be determined with greater certainty. This was called 
the design point, and is defined as the point toward lower 
solidity where the loading equals 1.10 times the minimum value. 
There is further justification in this step from the fact that 
the basic analysis neglects boundary layer. The blockage 
effect of a blade plus boundary layer must be greater than the 
effect due to the blade alone. Therefore, the true optimum 
solidity must be somewhat less than the apparent optimum based 
on non-viscous flow. The exact relation cannot be found with- 
out a detailed analysis of the boundary layer. This simple 
method sidesteps this formidable problem, and is at least 
qualitatively in the right direction. 

For the case of higher values of turning parameter and 
thick blades, there was some scatter in the data. Since the 
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data is all analytical, and should therefore be reasonably 
regular, some effort was made to find the source of the scat- 
ter. The relationship between the parameters cr" and C /S was 
suspected first. Comparative plots of these two parameters 
however, showed that they varied in a regular fashion with 
blade thickness. It was then decided that the scatter was due 
to the fact that velocities were measured only at 40 points 
along the airfoil. If the true maximum velocity occurred be- 
tween these finite points, then the value for boundary layer 
loading parameter would be in error. Remembering the fact 
that this error could be only on the low side was an aid in 
fairing the curves through the scattered points. 

It was noted that the data became more regular as the 
value of K approached one. For these airfoils, either the 
flow was more regular about the nose, with no sharp peaks in 
velocity, or the point of maximum velocity happened to coin- 
cide more often with one of the finite points. Examination 
of the pressure plots of Figure 4 points toward the former 
as the cause. With high values of K, the pressure changes are 
seen to be more gradual, and of lower magnitude. Because of 
this, the actual maximum velocity was always close to the 
velocity recorded at one of the 40 points. 

Figure 10 is an illustration of the effect of K in reduc- 
ing the boundary layer loading for a given value of turning 
parameter. With all other factors equal, the loading steadily 
decreased as the value of K increased toward unity. Note that 
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the effect of K was much greater xor the thick blades* For 
the case of extremely thin blades, K had no effect at all, 
which was not an unexpected result* 

4.3 Series 2 Plots 

Since the original series involved 30 graphs, some effort 
was made to summarize this data. The result of this effort 
was Figures 11 through 15, which collect all the data for a 
particular value of K on a single page. Boundary layer load- 
ing parameter and solidity are plotted along the coordinates 
as before, but now the solidity is design solidity. In effect, 
all the design points for a particular value of K have been 
plotted. Points of constant values of turning parameter have 
been connected, as have points of constant airfoil thickness. 
The result is that four interdependent parameters are shown 
for each point on the chart. Any chosen point then, repre- 
sents a particular value of blade thickness and turning para- 
meter, with the corresponding boundary layer loading and 
design solidity for the cascade. 

Figures 16 through 20 are plots made in exactly the same 
manner, except that optimum points were used instead of design 
points. It was necessary to change the scale for solidity 
along the horizontal axis, since optimum points represent 
higher values of solidity in every case. In addition some 
difficulty was experienced with irregularity in the data. 

This was due at least in part to the uncertainty in determin- 
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ing the value of solidity at the optimum point <> This optimum 
point data was included as a limiting case, to illustrate the 
maximum solidities to be used when the boundary layer is neg- 
ligibly thin. 

The choice of solidity as a coordinate is slightly mis- 
leading. It might be inferred from this choice that solidity 
should be used as a design input. However, solidity was 
plotted in this manner purely as a convenience in transferring 
data from the series one graphs. In design work, the two 
"inputs" could be boundary layer loading, and turning para- 
meter. The ,T output ,r information would then be blade thickness 
and solidity. An attempt was made to plot boundary layer 
loading parameter vs. turning parameter, but the resulting 
blade thickness and solidity lines were so nearly parallel 
that the information was difficult to read. 

4 .4 Comparison of NASA Separation Parameters 

The data on the series 2 plots contains all of the infor- 
mation necessary to compute both NASA separation parameters. 
This affords an opportunity to compare these two parameters 
directly. The boundary layer loading parameter S, and the 
NASA parameter have a constant relationship such that: 

Thus a Dj equal to 0.5 corresponds to a value of S of 4»25. 
This value is entered as a heavy line on all series 2 plots. 
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The approximation b, is related to the parameters B, B c , 
and c/S by the following relationship: 



Thus, since has been fixed at one value, each value of B 0 
and c/S will result in a particular value for D. This expres- 
sion for D shows that it cannot be a good approximation over 
a wide range of solidities, since it indicates that loading 
will always decrease whenever solidity is increased.. It has 
already been demonstrated that as solidity increases beyond 
the optimum point, loading increases also due to the blockage 
effect. 

Lines corresponding to a value of D=.5 and D=.6 have been 
plotted on Figures 11 through 20 as dashed lines. A note of 
caution is necessary in using the data on these figures. It 
must be remembered that there is really only one solidity 
shown on each chart, and that is the design (or optimum) soli- 
dity. The fact that solidity is shown as a variable is due 
solely to the fact that each combination of blade thickness 
and turning parameter demands a particular value of solidity 
at the design (or optimum) point. 

If the parameter D is a true approximation for D^, then 
the lines corresponding to D-.5 and D^~. 5 should be nearly 
the same. Examination of the figures shows that there is a 
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fair match for the case of optimum solidity* Generally, D 
give3 a value of loading which is Higher than D^, but the 
error i3 only around 10$, which is not unreasonable. For the 
case of design solidity, there is a greater discrepancy. For 
very thick blades, D still indicates a value of loading higher 
than D^, but this relation reverses so that for thin blades 
indicates the higher loading. 

The data in Reference 1 states that separation is likely 
to occur for values of greater than .5 or for values of D 
greater than .6. This in itself is an indication of a mis- 
match between D and Dj. If the occurrence of separation i3 
taken as the matching criterion, then the line corresponding 
to should match the D^=.5 line. For the case of optimum 

solidity it may be seen that the matching is very poor, where 
D indicates a much higher loading than D-^. At design solidity 
the matching i3 better, but D still indicates values of load- 
ing which are around 1S$ higher than D^. Extrapolating this 
trend leads one to suspect that, for some value of solidity 
less than design solidity, the two parameters D=.6 and D^=.5 
would probably match quite closely. From this, it is not un- 
reasonable to expect that the NASA obtained its experimental 
data from cascades whose solidity was less than the design or 
optimum as defined in this report. A few checks of source 
data for Reference 1 have verified this premise. Solidities 
were generally 1.0 for cascades with blade thicknesses of 
eight and ten percent. A check of the series 2 plots will 
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reveal that this is less than the assign or optimum point for 
any fairing shape. 

4»5 Comparison of Blade Thickness and Fairing Shape Effects 

Figures 11 through 20 clearly show the gains that can be 
made in turning parameter, at no increase in loading, by 
using thinner blades at the corresponding higher solidities. 
However, it was desired to compare the effects of changes in 
fairing shape and blade thickness on the same chart. Figure 
21 was made using the design point data, and Figure 22 was 
made with data from the optimum points. The value of S = 4*25 
was selected as the best estimate of maximum permissible blade 
loading. The values plotted were turning parameter vs. solid- 
ity at this constant value of blade loading. This resulted in 
lines which represented varying blade thickness at a constant 
value of K. Lines of constant blade thickness were then drawn 
in. It is interesting to note that these lines are nearly 
vertical. This means that the value of design or optimum 
solidity for a particular blade thickness is relatively inde- 
pendent of blade shape. This result was not directly predic- 
table from theory, and may prove useful in planning experi- 
mental investigations. 

As blade thickness is decreased, it is clearly seen that 
the importance of fairing shape decreases correspondingly. 

When blade thickness reaches about four percent, there is 
almost no discernible difference due to fairing shape. To 
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compare the relative effects of lairing shape and thickness 
changes, consider the following example., Select a blade thick 
ness of 12 # with K » .010, which is a relatively poor fairing 
shape with a blunt nose. For this shape at the proper solid- 
ity, a turning parameter of .74 is possible without separation 
If K * .990 the shape is nearly the best possible. 

Making this change increases the permissible turning 
parameter to .S3, or a 12# improvement. This same improvement 
could be obtained with the original shape by decreasing the 
thickness to a little less than 8# and increasing the solidity 
slightly. If the blades could be decreased to 4$, with the 
corresponding increase in solidity, the turning parameter 
could be increased to .92 without changing the loading. This 
is a 24 # improvement, or double the best improvement that 
could be obtained by fairing shape changes alone. As a limit- 
ing value, if the blades could be made infinite smally thin, 
the increase in turning parameter could be as high as 50#. 

It must be remembered that these performance figures are 
for cascades operating only at the ideal angle of attack. It 
may well be that these increases in peak performance are at- 
tained with some sacrifices in flexibility. This would be of 
no consequence in a constant output device , but for general 
application, the problem of operating off-design must be con- 
sidered. The analytical method described in this study can 
be modified to consider the case of operating at angles of at- 
tack other than the ideal. This matter is discussed more 
fully in the section on suggestions for further work. 
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5.0 Conclusions 

It is concluded that the method of generating cascades 
as described in this report is a workable one, and provides 
useful results. Cascades are generated as complete entities, 
which eliminates the need for separating the effects of the 
various design parameters. Use of a high speed digital com- 
puter means that large numbers of cascades can be generated 
and analyzed quickly. This allows exploration over a wide 
range of design parameters. 

Solidity in experimental cascades has often been arbi- 
trary. This investigation has revealed that the matter of 
solidity is of paramount importance. For any combination of 
fairing shape, thickness and camber, there is a particular 
value of solidity which will result in the lowest loading of 
the boundary layer. This loading can increase quite rapidly 
as the solidity is varied from the proper value. 

Data obtained in this investigation has yielded the un- 
expected result that the optimum value of solidity referred 
to above is relatively independent of the fairing shape of the 
blades. An important application of this result will be in 
planning experimental investigations of cascade performance. 

The actual values of solidity given in this report may be used, 
with the probability that they will be close to the optimum 
values for any sort of physical blade shape. 

Three parameters intended as indicators of the tendency 
toward flow separation were compared in this report. As a 
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result of this comparison, it is concluded that the parameter 
S, as defined in this report, is a useful indicator in this 
respect. It is further concluded that the value S « 4.25 is 
an acceptable design limit, which is as consistent with NASA 
data as this data is consistent with itself. 

With respect to the problem of increasing the performance 
of airfoils in cascades, it is concluded that the matter of 
blade thickness is of much greater importance than the fairing 
shape. The example was given of a cascade with blades of 12 $ 
thickness. Possible increases in turning parameter for this 
cascade due to improvements in fairing shape were shown to be 
on the order of 12$. Increases in turning parameter produced 
by decreasing the blade thickness to 4% with a corresponding 
increase in solidity, were shown to be 24$, or twice the pre- 
vious amount. The absolute limit of performance increase by 
this method would be given by the case of negligibly thin 
blades. For this case, the total increase in turning parameter 
would be 50$. 

It was pointed out in this study that increasing the 
performance of cascades by using thinner blades and higher 
solidities would probably reduce the flexibility of performance. 
The effects of operating off -design were not a part of this 
study since the method of generating the cascades produced 
airfoils operating only at the ideal angle of attack. It was 
stated that modifications could be made in the method to con- 
sider the off -design case, but this was not done due to time 
limitations. 
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6.0 Suggestions for Future Work 



This investigation represents only a beginning in the 
field of theoretical design. With the basic method and com- 
puter programs available, it should be a simple matter to ex- 
tend the scope of investigations similar to this one. Inves- 
tigations ought to be conducted for a number of average flow 
angles, so as to cover the whole range of possibilities in 
cascade design. Reference 2 for example, has some work on an 
average flow angle of zero, which represents the special case 
of an impulse turbine. There is no reason to prevent con- 
sidering negative average flow angles either. This would be 
the case of turbines rather than compressors, but the same 
method and computer program would apply. 

One might suppose that since the flow is accelerated 
through a turbine, the boundary layer experiences a negative 
pressure gradient and separation is not a problem. This is 
not exactly true, however. There are actually two pressure 
differences to consider in flow through a cascade. One is the 
effect perpendicular to the stream, which is a result of the 
lift being developed by each blade. The other is the longi- 
tudinal effect parallel to the stream which is the result of 
differences in pressure from the front to the rear of the 
blade row. Flow separation depends on the combination of these 
two effects. Thus, in decelerating cascades, the longitudinal 
effect reduces the allowable lift from each blade. In tur- 
bines, the longitudinal effect increases the allowable lift. 
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If the lift is increased enough to overcome the aid of the 
longitudinal effect, flow separation will still occur. A 
series of studies of this matter could establish the exact 
relation between these two gradients and the occurrence of 
separation. 

The analytical methods of this report would be a help in 
studying the occurrence of separation through solution of the 
boundary layer equations. The detailed pressure distributions 
provided by the computer at least open an avenue toward solu- 
tion of the boundary layer equations by direct integration. 

In any future work, it would be convenient to consider 
the velocity at more points near the leading edge of the blades. 
This would allow more accurate determination of the actual 
maximum velocity, and should eliminate some of the scatter 
found in the data for this investigation. The necessary modi- 
fications to the computer program are not difficult, and only 
lack of time prevented their incorporation in this investiga- 
tion. 

Another area of investigation would be a change in the 
basic transformation which would allow a finite angle at the 
trailing edge of the blades. This would make the transforma- 
tions similar to the Karmann-Trefftz transformations then, 
rather than the Joukowski. Eliminating the cusp or radius at 
the trailing edge then, would perhaps yield blade shapes which 
would be more physically useful. As pointed out earlier, some 
of the blade shapes developed in this investigation, while 
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mathematically useful, were physically poor because of the 
rounded trailing edge. These improved transformations would 
be particularly useful in the case of symmetrical impulse 
blading, where a finite angle is desired at both the leading 
and the trailing edge. 

As mentioned previously, it would be interesting and per- 
haps useful to consider the off -design case, where the blades 
are operating at some angle of attack other than the ideal. 

In principle, the method would be to generate the airfoils, 
then go back to the original complex potential and change the 
average flow angle. Keeping the same transformations would 
preserve the blade shape, but the circulation would have to 
be adjusted to keep the Kutta condition. The idea seems simple 
enough, but no doubt there will be many difficulties with the 
details. 

Still further work could be done in investigating other 
shapes for the camber line of the blades. Generally, the 
camber line for the blades developed in this investigation was 
the slightly flattened arc of a circle. The generation of 
other camber lines involves the expression for the flow poten- 
tial in the near circle plane. In this plane, a single doub- 
let generates the flow contour which transforms into one of 
the blades of the cascade. If instead, a series of doublets 
were used; by varying the strength and distribution of these 
doublets, great control could be exercised over the camber 
line and blade shape of the airfoils. The ultimate in this 
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direction would be where the input is a desired pressure 
distribution, with the proper camber and blade shape being 
automatically produced. The resultant family of cascades 
would be as important a milestone to cascade design as the 
NASA 6-series airfoils were to isolated airfoil design. 
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Boundary Layer Loading iarameter vs. Solidity 

Bor Constant Thickness Airfoils 
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Boundary Layer Loading Parameter vs. Solidity 

For Constant Thickness Airfoils 
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Boundary Layer Loading Parameter vs. Solidity 
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Boundary Layer Loading Parameter vs. Solidity 
Lor Constant i'hickness Airfoils 
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Figure 10. 

.Boundary Layer Loading Parameter vs. Solidity 
Showing the Effect of Parameter K Variations 
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Figure 11. 

Boundary Layer Loading Parameter vs. .Design Solidity 
For Various Cascades Beta = 45° K = .010 
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Boundary Layer Loading Parameter vs.Lesign Solidity 
For Various Cascades Beta = 45° K = .010 
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T’igure 12. 

Boundary nayer Loading Parameter vs. Lesign Solidity 
.For Various Cascades Beta = 45° K = .255 
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Boundary Layer Loading Parameter vs. .Design Solidity 
Bor Various Cascades Beta = 45° i£ = .500 
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Boundary Later Loading Parameter vs. Design Solidity 
for Various Cascades Beta = 45 iv = .745 
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Boundary Layer Loading Parameter vs. Design Solidity 
Bor Various Cascades Beta = 45° K = .990 
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Boundary Layer Loading Parameter vs. optimum Solidity 
Bor Various Cascades Beta = 45° K = .010 
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Boundary Layer Loading Parameter vs. Optimum Solidity 
i'or Various Cascades Beta = 45° & = .245 
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figure 20. 

Boundary Layer Loading Parameter vs. Optimum Solidity 
Por Various Cascades Beta = 45° K = .990 



54 



*4 3 % 2 M f? 







figure 21. 

Turning rarameter vs. Design solidity 

ior boundary nayer noading parameter of 4.25 

( Corresponds to LaSa Local Diffusion factor of .5 
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.figure 22. 

l'urning Parameter vs Optimum Lolidity 

for Boundary Layer Loading rarumetcr of 4.2b 

( Corresponds to JSJASa Local Lif fusion factor of .5 ) 
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APPPhuIX 



14ATHBMATICAL ANALYSIS UP Ti^ naTnuD uP jjaBlGh 

This appendix describes the general mathematical 
development of the method of designing airfoils in cascale. 
Some of the fine details, particularly on iterative proceea- 
ures, have been omitted, however, the entire calculation 
proceedure is given in the j^aivihAh program in the next appen- 
dix. This, along with the explanations given in the table 
of symbols, should allow analysis of the method to whatever 
extent as necessary. 

Basic relations in the hear Circle rlane 

The complex potential for the flow around the infinite 
row of closed bodies in the near circle plane may be written 
in the form; 

w = <f‘ + if = z + (a + iB) P-^(z) Ala 

with the corresponding complex velocity; 

dw = v - iv = e~^ + (A + iB) P~(z) Alb 

dz X y 

where 

P 1 (z) = K^x.y) + iJ 1 (x,y) = ± coth | A2a 

P 2 (z) = R 2 (x,y) + iJ 2 (x,y) = P 1 (z) 

= dz ^ coth f j = - t csch 2 | A2b 

By expanding these functions and separating tnem into 
their real and imaginary parts, it is found tnat: 



b7 



^(x^) 


, /sinh x | 


A3a 


2 Icosh x - cos y J 


Jj/x.y) 


= ' * ( 


sin y I 


A3b 


cosh x - cos y / 


* 2 U,y) 


1 

= * ? 


(cosn x cos y - l) 
(cosh x - cos y^ 


A3c 


J 2 (x,y) 




(sinh x sin y ) 


A3d 


= + T 


(cosh x - cos y)2 



xhe expressions for velocity potential, stream function, 
and velocity components now reduce to: 



$ = x cos ft + y sin ft 


+ AK x (x,y) - BJ-^Xjy) 


A4a 


y = - x sin ft + y cos ft 


+ AJ^(x,y ) + BR 1 (x,y) 


A4b 


V = cos ft 

A 


+ AU 2 (x,y) - BJ 2 (x,y) 


A4c 


Vy = sin ft 


- AJ 2 (x,y) - BR 2 (x,y) 


A4d 



i'hese relations represent the flow without circulation 
about an infinite series of nearly circular bodies spaced 
at equal intervals S = 2 Tf along the y axis, xhe contour of 
a typical body is shown in figure c j>. 

Location of Stagnation Points 

hue to the interference caused by thu presence of the 
infinite stack of bodies, the stagnation points on any body 
are displaced slightly from the position they would have 
occupied if the body were isolated in the flow, for this 
reason, some computation is necessary to locate these points. 



designating the front ana rear stagnation points by the 
subscripts a and d respectively, anu noting the type of polar 
symmetry involved, the following facts are observed. At the 
two stagnation points, both velocity components vanish by 
definition, ^lso, by reason of the odd symmetry involved, 

> y A = _ y B » and • However, since 

the points and B lie on the same streamline, ^ . 

a he only way in which these two conditions on f can be satis- 
fied is if f = - V'jp 0 . oince the stream function vanish- 

es at the stagnation points, it also vanishes at all other 
points on the body, since tney are on the same streamline. 

The velocity potential <f> on the other hand varies from 
a maximum negative value at A to a max: mum positive value at 
i3. 'i‘he value f B then defines the amplitude of the func- 
ti n on the contour of the body, ihis amplitude may be 
arbitrarily prescribed, and will then fix the size of the 
near circle contour in relation to the fixed spacing J = 2]f. 

Let y ) = and so forth, i’hen, at the point 

H from tne foregoing discussion: 



g = 

’ B 


cos + y B si np 


+ Aii lB - bo 1B - 


arbitrary value 
A5a 


V = 


-Xg sin^ + y B cos f 


+ A J iB + bh 1b 


= 0 


A5b 


v x 


cosf 


5 + ” 3J 2B 


= 0 


A5c 


v y = 


sin 1 


? - AJ 2 b ~ Bxv 2B 


= 0 


A5d 
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l^-st two oi the above equations may be solved simul- 
taneously for the unknown constants a and B with the results; 



- °os/2 + sin (3 

__ ~2 

n 2B + J 2B 

+ sin^3 + cos j2 

R 2 B 2 + J 2B 2 



A6a 



A6b 



substituting these results into the expressions for 4^ and V 
yields two equations in wnich the two unknowns are the 
coordinates and y^ of tne aft stagnation point. The 
equations are trancendental , however, so that an iteration 
technique is necessary for solution. Care was necessary in 
this area as all the functions are multi-valued in y, and the 
result could converge to a point differing by some integral 
multiple of 2 IT in the y dimension from the intended point. 

The general method involves starting with low values of 
in effect close to the isolated case where the location of 
the stagnation points :s known, uiiu working up to the desired 
value by modest increments. 



nocation of Points on the Contour 

^nce the location of tne aft stagnation point has been 
found, the constants a ana n can be evaluated from equations 
h6a and A6b . rhe other points on the contour can then be 
located from equations A5a and Abb by setting <p = 0 and 
assigning a suitable value to rp . 'i'nis computer program 
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divides the interval from -f B to s- into 20 equal parts, 
fhis yields 20 points along each surface of the airfoil, or 
40 points in all. As pointed out in the text, it would be 
helpful to subdivide some of these intervals near the leading 
edge to get a better check on the maximum velocity in this 
critical region. 

ihe equations for the location of these contour points 
are again trancendental, and a similar iterative technique 
is necessary for solution, wnce the coordinates of a point 
are known , the velocity components are found from equations 
a5c and .u5d. ihe final velocity v and the slope X of the 
contour at that point are found from the relations: 




v 

oC = arc tan A8b 

x 

At the stagnation points, this relation becomes indeterminate 
and L’Hopital’s rule is used for a solution. 

.Basic relations in the Circle rlane 

It is well knotm that the flow around an infinite row 
of arbitrary but identical contours can in principle be 
mapped into a corresponding flow about a single unit circle. 
If the flow is of zero circulation, and the direction of flow 
far from the contours is (3 then the corresponding flow in 
the circle plane is expressed by tnc- following complex po- 



tential and velocity: 
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du 

d $ 




e 



-i & 




splitting these into real and imaginary parts, and placing r 
=a = 1 leads to the result that along the unit circle, 



= 0 AlOa 

AlOb 

^ = constant 



ibis verifies that the circle is in fact a streamline, 
further, along this circle: 





cosh <£_ ± cos 
cosh/ - - cos & 



i 



+ sin 






2 arc tan 



( 



sin & j 

sinh J 



Alla 




cos/^cosh^ sin<? + s in ft 
cosh - 



sinh 6* cos 9_ 
cos 2 & 



Allb 



Singularities and stagnation points in the Circle ilane 
io find the stagnation point S, set = 0 . ihen 

equation Allb gives: 

tan <&£ = tan tanh a12 

Substituting this into equation aIIa produces a result which 
can be reduced to the form: 
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4 a = cos ^ in ^ — j + S±nf>.2 arctan ( ) nl3a 



wnere 



T = y^( cosh 2^ + cos 2ft) 



A13b 



ihese equations then must be solved for T then for 

assigned values of and which must be the same as the 

values for /3 and <^j^L previously employed in the near circle 

plane. Since the functions are trancendental, an iterative 

technique is again required. Newton's metnod is employed, 

J <Qo 

which involves calculation of the derivative z — — . he- 

# J r 

noting this by the symbol then: 



<?: = 



2 

-8'f 



B 






2T - cos 2i 



f) 2 - 



A14 



Applying the equations in this form however, requires an 
excessive number of significant figures in the evaluation of 
T . lo circumvent this problem, the following subltitution 
was used to eliminate T. 



1 = cos 






A15 



i'he x-'^evious equation then become 



Jf In/ 2 C0S 


/3 


+ A T 1 


— cos p -Ln / ^ 


i' 


J 



-] + Bia P - 2 arotan [f^r + ^ jl 



Al6a 
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. 







Al6b 




r 2 cos P 



+^ r i')^'x7 



the problem of finding T for assigned values of ft and*^ 
by newton's method has no difficulty with regard to signifi- 
cant figures. 1'he parameter^ 4 and its liyperbolic functions 
may now be found from: 



sinh^ 1 = y T ( 2 cos/3 + ( cosh 2^-1) 



A17 



'i'he initial estimate for hewton's method is obtained by 
neglecting 1 in the last term of equation al6a . the 
resulting equation may be solved explicitly for A 1 ♦ 



where 



A T = 



2 cos @ 

(J.i) 



I 




2 ft_ sin 
cos fi 



) 



AlOa 



A18 b 



these expressions provide an excellent first estimate for 
T, and in fact for the special case of^= 0 , they are 
exact. 

with ^3 and known, the location of any point ^ 
having an assigned value of is defined implicitly by 
e j.uation «lla. to establish the location of points corres- 
ponding to those in the near circle plane, all that is neces- 
sary is to have the respective values of bo the same in 
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the two planes, as before, an iterative solution is re- 
quired . 



Circulation 

'ihe complex potential ana velocity for a purely cir- 
culatory flow of unit circulation in the circle plane are 
as follows: 




dlAw) _ - 2i cos ft sinh f £ 

i 4 - 2 cosh/ 3 . J 2 +l 



A19b 



5 ■ 



= re 



i<s? 



A19c 



Separating these into real and imaginary parts, it is again 
found that the unit circle r = 1 is a streamline, along 
which the velocity distribution is given by 

AM - - sinh 2 Yj a 20 

cosh - cos 2 



This circulatory component must be evaluated at the same 
points # previously establishecykround the periphery of the 
circle. By superposition, the total velocity at these points 
is then given by: 




o “ 
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utagnaticn Joints with circulation in (Jircle xlane 



Just as the stagnation points it ana B are defined by 
the condition V = 0 , so are the stagnation points 1 and $ 
defined the corresponding condition U- 0 . introducing 

S' 

this condition into the foregoing equations yields the 
essential result 

e, &. = r = arcs in • sinh 2/ cos 

B 1 L A L 4 ■ c'osh^ cos T, 

A22 

once the points L and 1 have been located in the circle 
plane, the corresponding values of and may tie found 
from equation Alla, introducing these \alues into the near 
circle plane leads to the location of the points i» and 1 in 
that plane by solving equations A4a and A4b. 

Mapping between Circle and near Circle Hanes 

Mapping between the circle or $ plane and the near 
circle or z plane, as implied in tne preceeding discussion, 
is accomplished despite the fact that no explicit function 
for z in terms of § or vice-versa is known, -“-he mapping is 
done implicitly through the fact that the non-circulatory 
flows in the two planes must correspond, or in other words, 
through the condition that 

w\, z) = w { § ) A2pa 
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differentiating this gives 



dw 

dz 




A23b 



.But 



dw | = v 

dz 



■“•23 c 



and 




A23d 



Hence the mapping derivative becomes 



dz i = Vg» 
d^ ( v 



A23e 



xhe flew with circulation has been comxdetely defined 
in the foregoing development for the circle plane. In the 
near circle plane, however, only the non-circulatory flow 
has been explicitly defined, *.a explained in the text, add- 
ing circulation terms to the complex potential in the near 
circle plane produces unwanted distortions. Therefore, the 
method used is to map the known flow with circulation in the 
circle plane back into the near circle plane using the map- 
ping relations previously established. These mapping rela- 
tions are of course independent of circulation. thus, for 
the final flow with circulation 



w c U) = w c ($) = w {$) + B q Aw (5) 



A24a 
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differentiating this gives 




substituting in terms of velocities, arid introducing proper 
algebraic signs gives finally for the velocity in the near 
circle plane: 




V* -t- B^AV^a 



v 



= v + B 

o 




A25 



Basic delations in the Cascade Plane 

The basic idea of this entire method is to map the near 

1 

circles into a cascade of airfoils in much the same manner 
that the Joukowski transformation may be used to map a single 
circle into an airfoil. The equations for the transformation 
are identical in form to those already developed for the 
flow in the near circle plane, however, all coordinates are 
evaluated with respect to axes shifted slightly from the 
original position. Denoting these shifted coordinates by 
primes, the mapping function employed is specifically: 

J = e _l/ z' + (A' + iB') J^U') A25a 

-4^7 = U - iu = e" 1 ^' + (A' + iB') F 9 (z') A25b 

dz ' x y 2 
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Where 



A 26 a 



i'-^z • ) = coth 
i^U*) = -i csch 2 I 1 



A26b 



Hote that the above functions are exactly of the same form 
as equations Al and A2 . All subsequent results may readily 
be obtained through this analogy. 



determination of Airfoil chape 

In order to transform the near circle into the required 
airfoil shape, it is necessary to locate the singular points 
t and of the mapping transformation in the proper relation 
to the stagnation points n and i of the flow, normally, 1 
is located a short distance inside the contour, just opposite 
L on the surface. Similarly, point lies inside the con- 
tour on a line normal to the contour at T. nefer to i’igure 

5. J-he significance of these locations may be qualitatively 

PL 

understood from the fact that as the relative distance — 

= £ is made small, the radius of curvature at the nose of 
the airfoil decreases rapidly, and in the limit, as p = 0 
the nose becomes cusped. Similarly, the relative distance 
^ as pK controls trie tail curvature. In the limit as 

k b 

= 0 the tail becomes cusped in the manner of a eoukow- 
ski airfoil. l‘he origin of coordinates for the transforma- 
tion is shifted from the original location, the midpoint 
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between points A and B , to the new location, the midpoint 
between points P and Q. 



Velocity and iJlope p-long the Airfoil 

The complex potential including circulation in the cas- 
cade and near circle planes may be equated 

w (i*) = w (z ) A27a 

c * c 



differentiating 





A27b 



net 



( dv o \ = ve- 1 ' s 

Vd5"*J 


conjugate of velocity 
on airfoil 


A27c 


(!i). 


conjugate of velocity 
on near circle 


A27d 


(J^j- 




A27e 






*.27f 



Prom which, alter separating amplitudes and phase angles • 



V = 



2 / 



v + * 0 * v e(h) 

u 



A28a 




A28b 
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APrbNwiA b 



Fortran Program "Cascade" and Operating lnstructioxis 

on the following pages is tne FoxtldiAh program "Cascade" 
used in deriving tne data used in this study, xnput data is 
by means of data cards inserted at the end of the program. 

.my number of data cards may be used, x-he first card indica- 
tes how many additional data cards are to be read. 

nach additional card lists the number of values of the 
parameters , , x3 q , and x£ to be computed, their 

initial values, and the increment by which each value increa- 
ses. also listed are the print control quantities burn and 
NbhPh , which control the extent of the print out as indica- 
ted in the table of symbols, -the computer then calculates 
results for every possible combination of the five parameters, 
prints the results, and then proceeds to the next card. 

Ihe standard print out contains the five inpet parameters, 
the solidity and stagger angle of the cascade, the airfoil 
thickness, and the coordinates, slope, velocity, and boundary 
layer parameter for each of the 42 points of the airfoil. 

Condensed Print out 

when surveying large numbers of cascades, as was done 
for this study, it is convenient to condense the data to 
reduce the volume of output. Page 83 shows the modifications 
made to the basic program to accomplish this, briefly, all 
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of the standard print commands w ere removed, xhen commands 
to print the heading only were inserted early in the program 
right after changes in sigma. The modified print out then 
contains no airfoil data at all except for the maximum value 
of the boundary layer parameter, operation of the program 
is identical with the basic program. 

Iterative Scheme for Constant thickness Airfoils 

Page 84 shows further modifications made to the basic 
program to produce cascades with the desired airfoil tnick- 
ness. in effect this eliminates the parameter <£ since it 
is varied as necessary to produce the required thicknesses. 
The initial value of £ is still read from the data card, but 
the column which indicates the number of £ to be calculated 
no longer has significance. The program as written will 
demand thicknesses from 4 c /» to 24 °/<> in 45* increments, uther 
values for thickness can be obtained by modifying the appro- 
priate commands. No means of doing this by data cards is 
provided. The print out is of the condensed form described 
previously . 
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..JOB JENISTA 

PROGRAM CASCAOE 

OIMEMS^N XI(i.2{,YnU2).R I (l*2) i ,V!U2).A L I(42),AHI(42) t 0ELI(U2). T EI 

201 READ 204.18 
204 FORMAT (112) 

202 REAO 205. NBE 



rAI^TSFOR 2 ftiAR > CIRCLE 



BE J 



1AKI , nqpr -NixPR* NBZ,NEf> * NAK * BE1 * SG1,BZ1,EP1,AKl,BEI * SG1 ' BZ, » EBI * 
205 FORMAT (Sf2. 1CF5. 

PART I A CONSTAN 
BEJ*BE 1 

BE* l5t§93/1 80. ) 

SNBE=SINF(BE ) 

CNBE*COSF ( BE ) 

SGK*SG1 
003 K* 1 »NSG 
PHEB*SGK*3. 141593 
C22-SN8E 
C23*CN8E 

IF(K-I) 401,401,402 
0RUP*1.5708 
0ALP*0. 

SIP*0. 

RUP*0. 

ALP*BE 

IF(SGK-SIP-.S) 404,404,403 
S I G=S IP*.5 
RU=RUP*DRUP* .5 
ALU*ALP+DALP*.5 
IF(NEXPR) 52C, 520, 518 
WRITE OUTPUT TAPE 2 “ 

FORMAT ( / 7H ORUP 

1 5HSIG *F10.4, 5X _ 

2 5X 5HALU *F10.4/ ) 

520. CALL BAKER ( S IG, XU, YU, RU ,ALU, AU.BU, R 1U, AJ 1U, R2U, AJ2U , .000001 , 

A2U, ITU »C 22, C 23, NEXPR > 

DRUP* ( RU-RUP )/ .5 
DALP*(ALU-ALP)/.5 
SIP*SIG 
RUP*RU 
ALP*ALU 

I F(NEXPR) 522,522,521 

WRITE OUTPUT TAPE 2, 5 19, ORUP.OALP.SI P, SI G, RUP, RU, ALP , ALU 
CONTINUE 
GO TO 402 

RB*RUP*DRUP*(SGK-SIP) 

AL8*ALP+DALP*(SGK-SIP) 

IF ( NEXPR) 525,525,523 

WRITE OUTPUT TAPE 2, 524, DRUP ,DALP, SI P, SGK, RUP, RB, ALP , ALB 
7H DRUP *F8. 4, 5X 6HDALP * F9.4 



401 



402 

403 



518 

519 



>18 

>,519, DRUP, DALP, SIP, SIG,RUP,RU,ALP . 
'F8.4, 5X 6HDALP * F9.4, / 6H SIP 
5HRUP *F 1 0.4 , 5X 4HRU*F11.4, 5X 



ALU 

«F 9.4, 5X 
5HALP -F10.4, 



1 



521 

522 

404 



523 

524 



FORMAT( / 



LP * F9.I*, /' 6H SIP *F 9.4, 5X 
5X 4HRB *F 1 1.4, 5X 5HALP *F10.4, 



1 5HSGK *F 10 .4 • 5X 5HRUP *F10.4, 

2 5X 5HALB = F 1 6 • 4 / ) 

525 CALL BAKER ( SGK, XB,YB,R8, ALB, A, B, RIB, AJ1 B,R2B,AJ2B, .000001 • 

1 A2, IT, C22.C23, NEXPR) 

DRUP* ( RB-RUP ) / ( SGK-S IP) 

DALP* ( ALB-ALP) /(SGK-SIP) 

S I P*SGK 

RUP*RB 

ALP*ALB 

I F ( NEXPR) 527,527,526 

526 WRITE OUTPUT TAPE 2, 524, DRUP, DALP, SIP, SGK, RUP, RB, ALP , ALB 

527 CONTINUE 
SNALB*YB/RB 
CNALB*XB/RB 

ALBO *(180. /3. 141593)* ALB 

XA*-XB 

YA*-YB 



Tj 



ALA-ALB+3. 141593 

R A*RB 

X I ( 20)*XA 

X I (40) *XB 

Y I { 20 ) *YA 

YI(40)*YB 

R I ( 20 ) *RA 

R I ( 4 0 1 *RB 

ALI ( 20 1 *ALA 

ALI (401-ALB'*’ 6* 28 3186 

VI (201*0. 

VI(40)*C. 

R3B— 2.*(RlB*R2e-AJlB*AJ2B) 

AJ3B*-2.* ( AJ1B*R2B*R1B*A J2B) 

DVXB*A*R3B-B«AJ3B 

0VYB*B*R3B+A«AJ3B 

DVZB*SQRTF(DVXB*0VXB*DVYB*0VYB1 

AHI ( 40 ) =AS I NF ( DVYB/DVZB) -1.5707965 

AHI (20) -AHI (40) *3. 141593 

A22=A*A*B*B 

A23*R1*RHAJ1*AJ1 

A24=R2*R2+AJ2*AJ2 

W2ZB*2 »*SQRTF ( A22*A23*A24 ) 

IF (N0PR1120, 120,121 

!2°i WR a2 E ?¥ TPUT TAPE 2 * 122 ’ BEJ » SKBE » CNBE » SGK » PHEB * A » B » x0 » YB » Re » AL6C 

1220F0RMAT (1H1///_46X 28HRE SULTS IN NEAR CIRCLE PLANE // 



1 

2 

3 

4 

5 

6 
7 



10X 6HBETA -F6.2, 5H DEG. 6X 10HSIN BETA * F7.4, 8X 
10HCOS.BETA * F7.4, 6X 7HSIGMA » F7.3 t / 



1 0X 8HPHI ( B ) = F6.4, 9X 3n« * oa ->n d 

12X 1 H I 9X 1HX 1 1 X iHY 11X 1HR 5X IJhLAMBDA, 

. - SHIFT 



6X 1 1HALPHA, DEG. 3X5H ERROR 5X 10HITS 
12X 1 HO 3F12.4, F12.2, 6X 6H0.C000 8X4H* 
6H *1 

121 CONTINUE 

PART IB POINTS ON NEAR CIRCLE 
DC 12 1*1,19 
IF(I-I) 13,13,14 

13 PHN=.9*PHE8 
AL*ALB*.45 
R*RB 

GC TO 15 

14 A I* I 

PHN=PHEB*( 1.-AI/10.) ■ 

‘AL+( 3.141 593+ALB-AL 1/ (21 .-AI ) 



HA * F14.4, 8X 3HB « F14.4, // 

DEG. oX 1 HV 
// 

E12.2, 16, 



15 



49 

50 



36 

37 



AL J _ 

X=R*COSF(AL) 

XALL S SSR( X^Y,VX , VY,V,C22 ,C23, A, B, PHN, PHEB, ITS, SNALB, CNALB , 1 • , 
1 ITSIDE,ERS,NEXPR) 

XI ( I 1 *x 
YI ( I )*Y 

R*SQRTF ( X*X*Y*Y 1 
R I ( I )*R 
V I ( I 1 *V 

AHI ( I )*ATANF ( VY/VX) 

IF(VX)49, 50,50 

AHI ( I )*AHI ( I 1+3. 141593 

AHD TINU i(180./3. 1415931* AHI ( I ) 

AL*ATANF(Y/X 1 
IF(X) 36, 37.37 
AL*AL+3. 141593 
CONTINUE 

ALD *( 180. /3. 1415931* AL 

AL I ( I 1 =AL 

XI ( I +20 )*- ( XI ( I 1 1 

YI( I *20 1 ■— ( Y I ( I 1 1 

VIII +20)*- (V Kill ,, 

ALI ( I*20)*ALI( 1 1*3.14159 3 
AHI ( 1*20) * AH 1(1) 



74 






. 



r 








• - • • 



. 

' ! ’ • 

• •» *v 

■ ' 
















124 

125 
12 



528 

529 

530 



-5 1 IlX* Y» R » ALO . V • AHOtERS , I 



16 



R I ( I ♦ 2 0 ) * Rid) 

IF (NOPR) 124,124,12 
WRITE OUTPUT TAPE 2, 1 
FORMAT (I13,3F12.4, FI 
CONTINUE 
PART I I A CONSTANTS FOR CIRCLE 
IF(NEXPR) 530,530,528 
WRITE OUTPUT TAPE 2, 529 
FORMAT (///) 

I T3*0 

PCWER « (PHEe-2.*BE»SN8E l/CNBE 
DET*(2.*CNBE)/(EXPF( POWER 1-1.1 
ARG 1 « 1.*(2.*CNBE)/DET 
ARG2 » SNBE/(CNBE+DET) . . 

PHE* CNBE»LOGF( ARG 1 )*SNBE«2 .*AT ANF ( ARG2 1 
ER3*PHE-PHEB 
RER3*ABSF(ER3)/PHEB 
IT3* IT3+1 
T=CNBE+0ET 
Q1*CNBE«-T 
Q2=l .♦Ql*OET 

PHP*(-2./DET)»(T/Q1)«(T/Q2) 

DETNU * DET-ER3 /PHP 
IF(NEXPR) 516,516,505 
WRITE OUTPUT TAPE 2,506, 

FORMAT! 6H IT3 -I4.5X5H0E 
5X7HDETNU =E1C.4I 



ITS.ITSIOE 



505 

506 



l 



IT3, DET, PHE, RER3, DETNU . 

T «§10.4,5X5HPHE «F1 0.6, 5X6HRER3 



«E 10.4, 



) 

( (CH2G-1. 1/2. 1 
( (CH2G+1. 1/2. 1 



51 6 I F ( • 000001-R ER3 1 17,19,19 

17 I F ( I T3-25 1 18.19,19 

18 I F ( DETNU) 60,61,61 

60 DET=0ET/2. 

GO TO 16 

61 DET*DETNU 

GO TO 16 . 

19 CH2G=1.*2.«Q1*0ET 

SH2G=SQRTF (CH2G*CH2G-1 
SHG * SQRTF 
CHG = SQRTF 
C2=CH2G 
C3= SH2G 
C4*SHG 

TEB=ATANF( (C4/C5)*(C22/:23) ) 

TE 8 D *( 180. /3. 141 59 3 )• TEB 

CNTEB=COSF(TEB) 

TEA=TEB+ 3. 141593 

T E I ( 48 ) =TEB^ 6.283186 
'VTB s O. 

VTA=0. ' 

DVB=-(C3)/(C2-CCSF(2.*TEB) ) 

DVA=DVB 
VT I ( 20 1 =0. 

VTI (401*0. 

DVI ( 20 1 =DVA 
OVI ( 40) *DVB 
A22*C23*C5»CCSF(TEB) 

A23*C22*C4«SINF(TEB) 

A24=C2-C0SF(2.«TEB) 

W2TB*4.»( A22+A231/A24 

iOAnuBiTF 0 nnTPUT 6 UpE2. 127,SHG,CHG,DET,PHP.TEBD,DVB.RER3,IT3 

12TOFORHA T «/ £ te ]5 3THCORRESPOgOING ReSULTS s XM G £lRCl.| ^JnE/^ THOeL t , 



10X 12HSINH GAMMA * F8.4, 
E 1 0.4 , 8X 11HPHI PRIME 
1 5 X 1HI ?X 11HTHETA, DEG 
5HERR0R 11 X 3HITS // 1 5X 
El 9.2 , 113) 

’i 8 AR^ir?OINTS ON CIRCLE 



5X 12HC0SH 

*1 1X°7HV THETA 10X 13H0ELTA V THETA 
I HO F17.2, 15X 6H0.0000 F20.4, 



10X 
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DO 20 1=1,19 
A I * I 

PHN=PHEB*( 1.-AI/10.) 

IF l 1-1)21 ,21 ,22 

21 TEOLD = TEB 

TE = TEB ♦ SCRTF ( • 2»PHE B/W2TB ) 

GO TO 23 

22 TEOLD = TE 

TENU 1 = TE-.1*PHEB/VTE 
„ IF (TEA-TENU 1) 53,53,54 

53 TE = (TE+TEA1/2. 

GO TO 23 

54 TE = TENU 1 

23 I T4 »0 

25 C6=S INF ( TE ) 

C7*C0SF (TE 1 
C8*COSF(2.*TE) 

A2*L0GF ( (C5+C7 ) / (C5-C7) ) 

A3=2.*ATANF(C6/C4) 

PH=C23*A2+C22*A3 

A2=-(C23*C5*C6)*C22*C4*C7 

VTE=4.*A2/(C2-C8) 

DPH=PH- PHN 
IT4=IT4*1 
A2=ABSF(DPH)/PHEB 
IF(NEXPR) 511,511,512 

WRITE OUTPUT TAPE 2, 513, I, IT4, TE , CPH, VTE , A2 
FORMAT ( 6H I =14, 5X5HIT4 =14, 6X4HTE »E10.4, 



24 

26 

51 

55 

56 

57 
111 



5X4HA2 
111 
26, 111, 

51,51 ,55 



111 



512 

513 . _ .. . . ... 

1 5X5HVTE =£10.4. 5X4HA2 *E10.4) 
511 I F C I T4-25 1 24, 11 1. 

IF ( .000001 - A2) 

TENU = TE-DPF/VTE 
IF (TENU - TEOLD 1 
TE = (TE+TE0LD1/2. 

GO TO 25 

IF (TEA-TENU) 56,56,57 
TE = ( TE+TEA )/2. 

GO TO 25 
TE = TENU 
GO TO 25 

TED =(180. /3. 141593)* TE 



6X5H0PH =£10.4, 



TEI ( I) = TE 
VTI ( I ) =VTE 
DVI( I ) =-C3/ ( C2-C8 ) 

TEI ( 1+20) =TEI( I )+3. 14159 3 
VTI ( I*20)*-(VTI ( I ) ) 

DVI ( 1+20) = OVI ( I ) 

IF (NOPR) 129,129,20 
WRITE OUTPUT TAPE 2, 1 30, I , TED, VTE 
FORMAT (116, F 17.2, F21.4, F20.4, 
CONTINUE 

PART 1 1 1 A CONSTANTS FOR AIRFOIL 
BZL=BZ1 
D029L= 1 ,NBZ 
BAM=ABSF(BZL) 

I F ( 8AM-.00 1 ) 65,85,86 
TAU=0. 

X T = X R 
Y T= Y 8 
X I ( 4 2 )=XT 
YI (42)=YT 
VT=0. 

AHT=AHI (40) 

AHL=AHI (20) 

AHI ( 41 ) =AHL 
AHI(42)=AHT 
TET=TEB+6. 282186 
PHT*PHEB 
XL*XA 
YL*YA 



129 

130 
20 



85 



,OVI( I),A2.IT4 
E 19.2, 113) 
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X I ( 4 1 ) = XL 
Y I ( 4 1 ) =YL 
VL =0 . 

TEL=TE A 
PHL=-PHEB 
GO TO 87 

86 CONTINUE 
C2=CH2G 
C3=SH2G 
CUsSHG 
C5=CHG 
C22=SNBE 
C23=CNBE 

TAU = ASINF( ( 821-/4. )*(S H 2G/CHG)*(CN T EB/CNBEn 

TET=TEB-TAU+6.283186 

TEL=TEA+TAU 

A2 = S INF (TET) 

A3=C0SF ( TET) 

A4=(C5+A3)/(C5-A3) 

A5= A2/C4 

PHT=C23«*L0GF (A4)+C22*2.»ATANF(A5) 

A2=SINF ( TEL ) 

A3=C0SF ( TEL ) 

A4=(C5+A3)/(C5-A3) 

A5=A2/C4 

PHL=C23»LOGF(A4)+C22*2.*ATANF(A5) 
ENT=30.+10.«PHT/PHEB 
1 = 39 

32 A 1 = I 

3Q IF(ENT-AI) 3C, 31,31 
GO TO 32 

31 A2 = ( TET-TE I ( I ) )/ (TEI (1 + 1 )-TEI( I ) ) 

RT*RI (I)+A2«(RI( I + 11-RK ID 
ALT= AL I ( I ) + A2* (AL I ( I + l )-ALI( I>) 

XT=RT*COSF( ALT) 

YT*RT*SINF(ALT) 

ENL=30.+10.*PHL/PHEB 

1=21 

33 A I * I 

IF(AI-ENL) 34, 35,35 

34 1*1+1 
„ G0T033 

35 A2= ( T EL-T E I ( I ) ) / ( TE I ( I +1 I-TEI ( I ) ) 

RL = RI(I )+A2* (RI( 1+1 )-RI( I)) 

ALL = ALI (I)+A2*(ALI (I + D-ALimi 
XL=RL*COSF( ALL) 

'YL=RL*SINF(ALL) 

CALL SSR(XT,YT,VX, VY, VT, C22 ,C23 , A, B, PHT ,PHEB 
1 ITSIDE,ERST,NEXPR) 

X I ( 42 ) = XT 

Y I ( 42)= YT 

AHT=ATANF(VY/VX) 

IF ( VX)63, 64, 64 

63 AHT=AHT*3. 141593 

64 AH 1(42) =AHT 

CALLSSR (XL ,YL, VX, VY. VL ,C 22 ,C23 , A, B, PHL, PHEB, 
1 ITSIDE,ERSL,NEXPR) 

XI (41 )=XL 
YI (41 ) = YL 
AHL=ATANF(VY/VX) 

IF( VX)65»66, 66 

65 AHL=AHL+3. 14 1593 

66 AHI ( 4 1 ) =AHL 

87 CONTINUE 

TNBE=SNBE/CNBE 
TN8E2=TNBE-BZL/2. . _ 

C V2= ( 1 .+TNBE*TNBE)/( 1 .+T NBE2*TNBE2 ) 

CV*SQRTF (CV2 ) 

EPH-EP1 



»ITS»SNALB»CNAL8,-1., 



ITS,SNALB,CNALB,-1., 



! 
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DO 40 M=1 , NEP 
AKN= AK 1 

DC 41 N= 1 »NAK > 

IF(EPM) 99,96,99 

98 IF(N-I) 99,99,41 

99 CONTINUE 
XP=XL-EPM*RB»SINF< AHU 
YP*YL+EPM*RB«CGSF(AHl) 
XC*XT-AKN*EPP*RB*SINF< AHT) 
YQ=YT+AKN*EPM«RB*COSF( AHT) 

XZ=( XP+XQ) /2 . 

YZ= ( YP+YQ ) /2 • 

X GP = XQ- XZ 
YQP*YC— YZ 
X = XCP 
Y=YCP 
C2=S INH(X) 

C3=C0SH ( X ) 

C6=S I NF ( Y) 

C7=C0SF(Y) 

R1=(C2/(C3-C7) )/2. 

A J 1 = ( C6/ ( C3-C7 ) )/(-2.) 

PX=C3*C7-1 . 

PY=C2*C 6 

PZ=PX*PX+PY*PY 

R2=PX/(PZ*(-2.)) 

AJ2=PY/(PZ*2. ) 

QZ=R2*R2+AJ2«AJ2 

PU*Y*QZ-R2*AJUAJ2*R1 

PV=X»QZ-AJ2»AJl-R2*Rl 

BEP=ATANF(PU/PV) 

C22= S INF ( BEP ) 

C23=COSF ( BEP ) 

QX=AJ2*C22-R2*C23 

0Y=R2*C22+AJ2*C23 

AP=QX/QZ 

BP=QY/QZ 

AP22=AP*AP+BP*BP 

AP23=R1»R1+AJ1*AJ1 

AP24»R2*R2+AJ2*AJ2 

WP2Z8=2.»SGRTF<AP22«AP23*AP24) 

PART 1 1 I B POINTS ON AIRFOIL 
1 = 1 

42 X=XI ( I)-XZ 
Y=Y I ( I )— YZ 
C2=S I NH ( X ) 

C3=CCSH ( X ) 

C6=S INF ( Y ) 

C7=C0SF(Y) 

Rl=C2/( 2.* (C3-C7 ) ) 
AJl=C6/(-2.*(C3-C7)) 

PX=C 3*C7- 1 • 

PY=C2*C6 
PZ=PX*PX+PY*PY 
R2=PX/(PZ*(-2. ) ) 

AJ2=PY/(2.*PZ) . . 

XI I ( I )=X*C23-»Y»C22 + AP*R1 — 8P« A Jl 
ETI ( I )=-X*C22+Y«C23+AP*AJWBP«Rl 
A2=C23+AP*R2-BP*AJ2 
A3=C22-AP*AJ2-BP*R2 
A4=SQRTF( A2*A2+A3*A3) 

I F ( A4— . 0001 ) 171,171, 172 
171 RP3=-2.*(R1»R2-AJ1*AJ2) 
AJP3=-2.*( AJ 1 *R2+R 1 *AJ2) 
DUXP=AP*RP3-BP*AJP3 
DUYP=BP*RP3+AP*AJP3 
DUP = SORTF ( DUXP*CUXP+DUYP *DUYP ) 
EPI ( I)=ASINF(DUYP/DUP)-1 .5707965 
IF(ALI( I)-ALI(30) ) 173,173tll3 
EPI (I)=EPI ( I >+3. 141593 
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173 



172 

600 

113 



88 

89 

38 

39 

43 

44 



59 

62 

46 

47 



710 

716 

715 



GO TO 113 

EPI< 1 1 = AS I NF (A3/A4) 

IF(A2) 600,1 13,1 13 
EPI ( I 1 = 3. 141 593-EPI ( I ) 

46=lBSF(VTIlin ,,4 ' 5,3, * UH, " , ' E,,,,,n 
IF(A6-.0001) 88,88,89 
A5=SQRTF(W2Ze/W2TB) 

GO TO 38 
A5«vim/vTim 
IF(A8SF(A4)-*0001 ) 39,39,59 
I F C 8AM-.0001 ) 43,43,44 
VVI ( I l-SQRTF < W2ZB/WP2ZB1 
GO TO 62 

R3=-2.*(Rl*R2-AJ1*AJ21 
AJ3=-2.*(AJ1 »R2+R1*AJ2) 
WP2Z=SQRTF(AP22*(R3*R3+A J3*AJ31 ) 

SHX = S INH( X I ( I ) ) 

CHX = COSH( XI ( D) 

SNY=SINF(YI { I)) 

CNY=COSF ( Y I { I) ) 

R 1 * ( SHX /2 • ) / ( CHX-CNY ) 

AJl=(-SNY/2. ) / (CHX-C NY) 

PX=CHX*CNY-1 . 

PY=SHX*SNY 

PZ=PX*PX+PY*PY 

R2=PX/(-2.*PZ) 

AJ2=PY/(2.*PZ) 

R3*-2.* (R1 *R2-AJ 1 * A J2 ) 

AJ3=-2.*( AJ1*R2+R1*AJ2) 

W2Z=SQRTF( (A*A+B*B)*(R3*R3+AJ3*AJ3) ) 

SNTE=SI NF { TE I ( I ) ) 

CNTE=COSF ( TE I ( I ) ) 

SN2TE=SINF(2.*TEI (I) 1 
CN2TE=C0SF(2.«TEI( I ) ) 

B22=CNBE*CHG*CNTE 

B23=SNBE»SHG«SNTE 

B24=CH2G-CN2TE 

W2T=-4.*( 822+823 )/B24 

DW2T*( 2. *SH2G*SN2TE)/< 624*8241 

A25=W2Z-A5*A5*W2T 

A26=(DVI( I 1 * A 25 1 /A6-A5*A 5*CW2T 

VVI ( I ) = ( W2Z+BZL* A26 1/WP2Z 

GO TO 62 

VVI ( 11= (VI (I )+BZL*OV 1(11 *A51/A4 
I F ( 1-42 1 46,47,47 
1 = 1 + 1 
GOT 042 
CONTINUE ' 

DELX=XI K421-XI I (411 

DELY=ETI (421-ETI (41 1 

CH0RD=SQRTF(0ELX*0ELX+DELY*0ELY1 

SOLI 0=CH0RD/6 .283186 

SNA=OELY /CHORD 

CNA=DELX/CHOPD 

AROT=( 180./3.141593)*ASINF (SNA1 
STAG=( 180./3 . 14 1593 1 * (BE P+AS I NF ( SNA1 1 
AHTO =( 18C./3. 1415931 * AHT 

TETD =( 18C./3. 1415931 • TET 

AHLO =(18C. /3. 1415931 • AHL 

TELD =(18C./3. 1415931* TEL 

TAUD ={ 180. /3. 1415931* TAU 

BEPD =( 180. /3. 1415931 • BEP 

I AA= 20 

DO 717 IDO= 1 , 1C 
IAA =1 AA-1 
1= 35 

TDD=T ANF ( DEL I ( I AA 1*3. 14 1 593/ 180. 1 

AA= XNI ( I l-XNI ( I AA1 - (Y NI( I AA1-YNI ( 1 1 1 *TOD 
IF (AA1 714,714,713 
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713 



71 4 



717 



718 



)• 



1 = 1 - 1 
GO TC 715 

BRE=ATANF UYNMI + 1) -YM I(I))/(XNI(I*1)-XNI( I ) ) 1 
DD6L = .5*(BeE+OELI( IAA)*3. 141593/180. ) 

CC = ABSFI IYNI ( IAA1-YNK 1*1 ) )/CCSF(OCEL> ) 

DDI(ICp) = CC+ (XNKI + D- XNI ( I A A > -IYNI(IAA) - YN I ( I ♦ 1 ) 

1 TANF(DDEL) )*SINF ( BBE ) 

WAX I E = ODI ( 1) 

DO 718 I DO = 2,10 
WAX I E = MAX 1 F (WAXIE, 00 I ( I D D ) ) 

WRITE OUTPUT JAPE 2,302, BEJ,SGK,BZL,EPM,AKN, STAG, SOUO. WAXIE 
302 FORMAT ( 1 Hi //52X 15HC ASCAOE RESULTS// 

18H DESIGN PARAMETERS/ 

9H, BETA *F6.2, 5H DEG. 6X7HSIGMA ■ F6.3,10X8HB ZERO ■ 

F6.2, 9X 9HEPSIL0N = F5. 3,9X7HKAPPA « F6.3,// 

17H OTHER PARAMETERS/ 

1 1H STAGGER*F6.2,5H DEG. 4X 10HS0LIDITY -F6.3 ,7X 17HMAX. 
6CKNESS = F5.3) 

IF (NOPR) 140,140,141 

140 WRITE OUTPUT TAPE 2 , 304, XB.YB, PHEB, XT, YT . AHTD, TETD, PHT , XL , YL, AHLD, 
1 TELD, PHL,XP,YP,XQ, YQ .T AUO , XZ, YZ, AP. BP.BEPD 
3040F0RMAT (3X4HXB =F8 .4 , 1 IX 4HYB *F0 • 4 , 1 1X8HPHI ( B ) »F8.4,/3X 



1 

2 

3 

4 

5 



THI 



1 

2 

3 

4 

5 

6 

7 

8 



ro«H f i lAunrnnoi •r u* *♦ f / . 

4HXT = F8.4.11X4HYT =F8.4, 1 1 XI OHALPHAI T ) * F6.2. 5H DEG. 
2X1 OHTHETA ( T ) = F6.2, 5H DEG. 2X8HPHI ( T ) * F8.4./3X 
4HXL = F8.4. 11X4HYL * F8.4, 1 1 X 1 OHALPHA ( L 1 « F6.2, 5H DEG. 
2X10HTHET AIL) = F6.2, 5H DEG. 2X 8HPHKL) * F8.4./3X 
4HXP =F8.4, 1 1X4HYP =F 8.4 , 1 1X4HXQ *F8 .4, 1 1 X4HY0 *F8.4,11X 



F8.4, 11X4HBP * F8.4.11X 



12X9HI V/V2)»»2//1 



1XP =F8.4. 1 1X4HYP =F 8.4 , 1 1X4HXQ 
5HTAU =F 1 3. 2 , 5H DEG.,/3X 
4HXZ = F8.4, 11X4HYZ * F8.4, 11X4HAP 
12HBETA PRIME = F6.2, 5H OEG. ) 

141 WRITE OUTPUT TAPE 2,303 

3030F0RMAT (//42X 13HAIRF0IL SHAPE 32X21HVEL0CITY DISTRIBUTION/ 

1 88X 19HAT DESIGN CONDITION/ 

2 . 1 0X 1 H 1 1 7X 1 HX 19X 1HY 1 4X1 1HSL0PE, DEG. 1 1 X6HI V/ V2 ) 

1 = 1 

DO 90 1*1,42 

PRX=( XI I ( I )-XI 1(41) 1/CHORO 
PRY*CETMI)-ETI(41 )) /CHORD 
XNI ( I )=PRX*CNA+PRY*SNA 
YNI ( I ) = PRY*CNA-PRX*SNA 
DELI ( I ) =OEL I (I 1-AROT 
VVI ( I ) =CV*VV I ( I ) 

V2I ( I) = VVI(I)*VVI( I) 

90 CONTINUE 

WRI TEOUTPUTT APE2 ,305 , ( I, XNI ( I ) , YNI ( I ) ,DEL I ( I ) * VVI ( I ) ,V2I( I), 

305^FORMAT^ { II 1, 2F20.4, F20.2, 2F20.4) 

PARTIV PARAMETER CHANGES 
48 AKN=AKN+AK I 
4>1 CONTINUE 
EPM=EPM+EBI 
40 CONTINUE 
BZL=BZL+BZI 
29 CONTINUE 
SGK= SGK+SG I 
3 CONTINUE 
BEJ=8EJ+BE I 
2 CONTINUE 
I D=I D+l 

IF(IB-ID) 202,203,202 

^PART^V ITERATION SUBROUTINE FOR POINTS ON NEAR CIRCLE 
PART v i.tKAMLJN # V “ C2 2 , C 2 3 , A, B, PHN, PHEB , I TS , SNAL B , CNALB , 



SUBROUTINE SSR (X«Y.VX, 
1 SIDE , ITSI DE, A2 , NEXPR ) 



517 

510 



1 



I TS = 0 
ITSI DE=0 „ 

I F(NEXPR) 8,8,517 

F0RMAT?//5H T ITS^JhX 9X1HY 9X2HR1 7X3HAJ1 8X2HR2 7X3HAJ28X 1HA 
9X1HB 8X2HVX 8X2HVY 8X2HA2 5X4HSI0E/) 
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508 

509 
507 

9 

10 



8 C 2= S INH ( X ) 

C3=C0SH { X ) 

C6=S INF (Y) 

C7=C0SF(Y) 

A2=C3-C7 
R1 *C2/( A2*2. ) 

AJ1=C6/ (A2*{-2. ) ) 

PX=C3*C7- 1 • 

PY=C2*C6 
PZ=PX*PX+PY*PY 
R2=PX/(-2.*PZ ) 

AJ2=PY/(2.*PZ) 

PH=X*C23+Y*C22+A«R1-B*AJ 1 

PS=-X*C22+Y*C23+A«AJ1+B*R1 

DPH* PH— PHN 

VX=C23+A*R2-E*AJ2 

VY=C22-A*AJ2-B*R2 

V=— SI0E*SQRTF( VX* VX+ VY*V Y) 

A2=(ABSF(DPH)+AESF(PS) ) / PHEB 
I TS= I TS+1 

IF(NEXPR) 507 , 507 *508 

£nJJ£T°yT£ UT 2,509, IJS,X.Y,R1,AJ1,R2,AJ2,A,B,VX,VY,A2,ITSIDE 
FORMAT (15, 10F10.4, E10.4, 15) 

IFUTS-25) 9,100,100 
IF ( .000001— A2 ) 10,100.100 
0X=( VX*OPH— VY*PS ) / (V*V I 
0 Y= ( VY* DPH+VX »PS ) / ( V*V ) 

X=X-DX 
Y=Y-DY 
GO TO 8 

YY=Y*CNALB-X«SNALB 
IF (YY-SIDE) 101,11,11 
ITSIDE=ITSIDE+1 

- 102,11,11 



100 
101 

I F ( ITSIDE-3) 

102 YY=-YY 

xx=x*cnalb+y*snalb 

X=XX*CNALB-YY*SNALB 
Y=YY*CNALB+XX*SNALB 
GO TO 8 
11 RETURN 
END 

PART VI ITERATION SUBROUTINE FOR INITIAL STAGNATION POINTS ON NEAR CIRCLl 
SUBROUTINE BAKER ( SGK ,XB , YB ,RB, ALB , A ,B ,R 1 , A J 1 ,R2 , A J2, ERM , A2 , 

1 I T,C22,C23,NEXPR) 

IF(NEXPR)515, 515,5 00 

500 WRITE OUTPUT TAPE 2,501 

501 FORMAT ( //5H IT 4X2HXB 8X2HYB 8X2HR1 7X3HAJ 1 8X2HR2 7X3HAJ2 8X1HA 
1 9X1HB 9X2HPS 8X2HDP 8X2HA2/) 

515 CONTINUE , 

I T= 0 

PHEB=SGK*3. 141593 
SNAL B= S INF ( ALB ) 

CNALB= COSF(ALB) 

XB=RB*CNALB 
YB=R8»SNALB 
4 C2=S INH ( XB ) 

C3=C0SH ( XB ) 

C6=SINF{ YB) 

C7=C0SF ( YB ) 

A2=C3-C 7 
R1=C2 / ( A2*2 • ) 

AJl=C6/(A2*(-2. ) ) 

PX=C3*C7-1 . 

PY=C2*C 6 
PZ=PX*PX+PY*PY 
R2=PX/( PZ*(-2. ) ) 

AJ2=PY/(PZ*2.) 

QX=AJ2*C22-R2*C23 

QY=AJ2*C23+R2*C22 

QZ*R2*R2+AJ2«AJ2 
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B=QY/QZ 

PH=XB*C23+YB*C22+A*R 1-B* AJl 

PS=-XB*C22+YB*C23+A«AJ1*B*R1 

0P=PH-PH6B 

A2=( ABSF( DP) ♦ABSF(PS) J/PHEB 
ITMT+1 

IF(NEXPR) 51 4.514,503 

503 WRITE OUTPUT TAPE 2 . 504, I T, XB, YB ,R1 , A J1 ,R2, A J2 , A ,B ,PS ,0P, A2 

504 FORMAT ( 1 5 , 1 CF 1 0 . 4 , 1 E 10. 4 ) 

514 CONTINUE 

I F { I T-25 ) 5.7,7 

5 I F ( ERM-A2 ) 6,7,7 

6 A2*C3-C7 

R1X=(1.-C3*C7)/U2*A2*2. ) 

RlY=(C2*C6)/(A2*A2*(-2.) ) 

A J 1 X =— R 1 Y 

A J 1 Y=R 1 X 

PXX=C2*C7 

PXY=-C3*C6 

PYX=-PXY 

PYY=PXX 

PZX=2.*( PX*PXX+PY*PYX) 

PZY=2.* (PX*PXY+PY*PYY) 

R2X={PZ*PXX-PX»PZX)/I (PZ*PZ)*(-2.)) 

R2Y= {PZ*PXY-PX«PZY)/(PZ*PZ*{-2. ) ) 

A J2X =-R 2Y 
AJ2Y=R2X 

QXX=-R2X*C234AJ2X*C22 
QXY=-R2Y*C23+AJ2Y*C22 
QYX=QXY 
Q YY=-QXX 

0ZX=2.* (R2«R2X+AJ2*AJ2X) 

QZY=2.*{R2«R2Y+AJ2*AJ2Y) 

AX=(QZ*QXX-QX*QZX)/(QZ»QZ) 

A Y= ( QZ* QXY-QX*QZ Y ) / ( GZ*QZ ) 

BX={QZ*QYX-QY*QZX)/(QZ*QZ) 

BY=(QZ*QYY-QY*CZY)/(QZ*QZ) 

PHX=C23+A*R1 X+AX*Rl-B*AJ 1X-BX*AJ1 
PHY=C22+A*R1Y+AY»R1-B*AJ 1 Y-BY*A J V 
PSX=-C22+A*AJ1X+AX*AJ1+B»R1X<-BX*R1 
PSY=C23+A*AJ 1Y+AY«AJ1+B*R1Y*BY*R1 
A2*PHX*PSY-PSX*PHY 
DX= { PSY*DP -PHY«PS)/A2 
DY=(PHX*PS-PSX*DP )/A2 
XB=XB-DX 
YB=Y 8-DY 
.GO TO 4 

7 RB=SQRTF(X8*XB*YB«YB) 

ALB=ATANF{ YB/XB) 

IF{XB)28»45,45 

28 AlB=ALB+3. 141593 
45 CONTINUE 
RETURN 
END 
END 



5 6 1 5 45. 2.5 



5 .001 .01 



• 4 



.15 0 .245 1 0 



.END 
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PART IA CONSTANTS FOR NEAR CIRCLE 

BE J*BE1 

002 J = 1 jNBE 

BE= (3.141593/180.) • BEJ 
SNBE=SINF ( BE) 

CNBE=C0SF ( BE ) 

SGK=SG1 

003 K-1«NSG 

740 WRITE OUTPUT TAPE 2# 741, BEJ, SGK 

741 FORMAT ( 1 Hi /26HCA SCAOE RESULTS FOR BETA »F6.2,5H DEG. 4X7HS IGMA » 

1 F6. 3, //75H B ZERO KAPPA EPSILON SOLIDITY THICKNE 

2SS BL PARAMETER/) 



2. 10 

(WAXIE, 



00 718 IDD 
718 WAXIE = MAX1F 
751 PRAM 5 * V2I (1) . 

DO 760 1= 1, 19 
760 PRAM * MAXI F (PRAM 



ODI( IDD)) 



PRAM * MAXI F (PRAM. V2I(I)) 

WRITE OUTPUT TAPE 2, 755 , BZL* AKN. EPM ( 
FORMAT ( 3X,F6.2*2F 10. 3*2F14.3,F 16.4) 

PARAMETER CHANGES 



755 . ... 

PARTIV 



SOLID, WAXIE, PRAM 



riO JJli? 1 4 . CJ41 i’i uw o iuiv COi«o'i\^..x 



07 CONTINUE 

TNBE=SNBE/CNBE 

TNBE2*TNBE-BZL/2. 

CV2»( 1.+TNBE+TNBE l/C 1.+TN6E2+TNBE2) 
CV*SQRTF(CV2) 

. AKN*AK1 
00 41 N=1 »NAK 
WAXIT * 0 
THIK * .04 
EPM*EP1 

IF(EPM) 99,98,99 

98 IF(N-I) 99,99,41 

99 CONTINUE 



718 

756 

757 
761 

764 

765 



763 



751 

760 



1 TANF(DDEL) )*SINF ( BBE I 
WAXIE = 001(1) 

DO 718 100 = 2, 10 

WAXIE = MAX IF (WAXIE, DDIUDD)) 

IF (WAXSG) 757, 757. 751 
IF (THIK - WAXIE) 761, 751, 763 
IF (WAXIT) 764, 764, 765 
EPM = EPM - .01 
GO TO 99 

EPM = EPI ? ” *01*1 (WAXIE - THIK ) /(WAXIE - WAXIT)) 

WAXSG = 100.0 

GO TO 99 

EPM * EPM ♦ .01 

WAXIT = WAXIE 

GO TO 99 

PR AM= V 2 1 ( 1 ) 

DO 760 1= 1, 19 
PRAM = MAXI F (PRAM, 



WRITE OUTPUT TAPE Zt 755,. BZL, AKN, EPM, SOLID, 



, V2I ( I ) ) 

••■•--- -Ulrui I«rc 2, 755, BtL, AIMV, 

755 FORMAT ( 3X, F6.2,2F10. 3,2F14.3,F16.4) 
EPM=EPM+EBI 
THIK = THIK ♦ .04 
WAXIT = 0 
WAXSG = 0 

IF (THIK - .25) 99, 99, 48 
PARTIV PARAMETER CHANGES 
48 AKN*AKN+AK I 
41 CONTINUE 
BZL*BZl+BZI 
29 CONTINUE 
SGK=SGK+SGI 
3 CONTINUE 
BE J=BEJ+B6I 
2 CONTINUE 
10 * 10+1 

, IF(IB-IO) 203,203,202 
203_ END 



WAXIE, 



PRAM 
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